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ABSTRACT 



The sea breeze circulation was investigated using a combination of acoustic doppler 
sodar, doppler lidar and conventional observations in central California on the coast of 
Monterey Bay in September 19S7. The study was called LASBEX (Land Sea Breeze 
Experiment) and used the combined effort and resources of the Naval Postgraduate 
School (XPS), NOAA Wave Propagation Lab (WPL) and Naval Environmental Pre- 
diction Research Facility (NEPRF). 

A monostatic three-axis phased-array doppler sodar was able to detect the various 
features of the sea breeze front, which has many of the characteristics of a gravity head 
current produced in laboratory experiments, such as the longitudinal vortex, lobe and 
cleft structure, and strong regions of turbulence and convergence. Profiles of the wind 
field and its component variances were generated even - five minutes. Vertical velocities 
up to 1.0 m's were detected in the updraft region at the front. Vertical velocity fields 
were used to derive the horizontal divergence associated with the sea breeze and values 
of convergence in the range of 4 x 1 0 _3 s _1 and greater were calculated. These are believed 
to be more accurate values than any found in previous sea breeze studies. 

Doppler lidar radial velocity information was useful in mapping the various layers 
and wind structure of the sea breeze. The wind field of the sea breeze seen by the lidar 
was compared with the sodar. The results of this comparison showed that the sodar was 
well suited to investigate features under 600 meters in height, such as the updraft region 
behind the front and the variance of the wind, whereas the lidar was better at describing^ 
the larger scale flow, such as the return flow of the sea breeze, the convergence zone 
along the sea breeze front and prevailing synoptic winds. Together, these instruments 
gave a comprehensive and complete description of the sea breeze frontal structure, evo- 
lution and movement. 
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INTRODUCTION 



A. GENERAL DESCRIPTION 

The land-sea breeze is a diurnal mesoscale circulation that occurs on sea coasts due 
to the difTerences in surface temperatures between the land and sea. A typical horizontal 
temperature gradient is l°C/20 km (Atkinson, 1981). This situation induces a thermally 
direct circulation that alternates between onshore or sea breeze and offshore or land 
breeze. The circulation occurs during the day due to solar heating of the land which 
warms faster that the sea. This results in a horizontal pressure gradient with slightly 
higher pressure over the land. The pressure differential causes air to flow at higher levels 
towards the sea, thereby leading to convergence aloft. The resulting counter-flow at the 
surface is the sea breeze. Fig. 1. shows a lake breeze circulation, which has a similar 
structure as the sea breeze circulation. 

The land breeze is due to the same thermal structure except it is reversed. At night, 
the land loses heat quicker than the sea and, therefore, sets up a circulation opposite to 
the sea breeze. In both cases, the flow at the surface is from the cooler region toward 
the warmer region. By comparison, the sea breeze is the stronger circulation, since the 
heating of the land provides for a more intense heat source to drive the flow. 

The sea breeze occurs more frequently given a weak horizontal pressure field 
(Estoque, 1962: Simpson et ah, 1977). This is because the thermal effects have to alter 
any existing conditions. If the existing pressure field is strong enough, the resulting 
gradient winds can overwelm the land-sea breeze effect. 

Under some conditions, an offshore gradient wind will generate a frontal zone at 
the leading edge of a sea breeze (Atkinson, 1981). Estoque (1962) notes in his model 
that an offshore flow counter to the sea breeze will enhance the horizontal temperature 
gradient. The sea breeze front is also characterized by gradients in humidity, wind speed 
and wind direction. This sea breeze front is similar in form to a synoptic scale cold front, 
because the cooler moist marine air is advancing into the warmer dry land regime. 
Passage of such a front is a measurable feature and can induce significant mesoscale 
meteorological phenomena such as showers, thunderstorms and gusty winds (Pielke, 
1974; John, 19S0; Richiardone and Pearson, 1983). The horizontal extent of the sea 
breeze is large enough that it can develop a component parallel to the coast due to 
coriolis deflection (Defant, 1951; Walsh, 1974). 
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Fig. 1. Typical Lake Breeze Circulation, (from Keen and Lyons, 1978). 
B. MOTIVATION FOR THE STUDY OF LAND-SEA BREEZE 



The land-sea breeze circulation and sea breeze front can affect the dispersion of 
aerosols and pollutants over a wide area. This effect has been observed and modelled 
many times (Lyons, 1972; Lyons and Olsson, 1973; Anthes, 197S; Rondo and Gambo. 
1979; Kitada et al., 19S6; Ogawa et al., 19S6). In a sea breeze, particulates and gases 
can be carried inland over long distances. Collis and Russell (1976) summarized some 
of the important effects that particle pollutants have on climate, such as cloud forma- 
tion, precipitation and transfer of solar radition. Pollutants can enhance cloud forma- 
tion and result in acidic rainfall. Particulates and gases can block incoming solar 
radiation and trap outgoing long-wave radiation. 

In the case of lidar (light detection and ranging) sensors, concentrations of 
pollutants such as X0 2 and S0 2 can attenuate or distort the return signal so that infor- 
mation may be misinterpreted (Measures, 19S4). 
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Measuring the intensity and mapping the structure of the sea breeze can give a 
better understanding of how particles are transported and dispersed. Satellite data could 
be used to detect a sea breeze occurrence. In addition, forecasts of pollution events may 
be feasible. Surface data also can be used to evaluate satellite observations of the 
planetary boundary layer and the mesoscale features that occur during the sea breeze. 

C. PREVIOUS STUDIES 

A study was conducted by Johnson and O'Brien (1973) of the sea breeze along the 
Oregon coast. In their study, the sea breeze was found to have a maximum depth of 700 
meters, with a maximum speed of 7.5 m s at the surface. The zonal wind was utilized 
to describe the variability of the sea breeze because of the north-south orientation of the 
Oregon coast. The inla; J penetration of the sea breeze front was up to 60 km. A return 
flow was found between 1000 and 1500 meters. Johnson and O'Brien made the following 
observations: the low-level onshore flow of the land-sea circulation was located entirely 
within the marine layer; the marine layer tends to deepen briefly at the onset of sea 
breeze due to the low level convergence ahead of the sea breeze front; the return flow 
above the inversion appears to be in the form of surges, probably in response to the 
surges observed in the sea breeze; and a distinct wind maximum follows the sea breeze 
front inland. 

Simpson, et al. (1977) observed sea-breeze fronts that penetrated over 45 km inland 
with some fronts moving as far as 80 km. These fronts had an average advance speed 
of 3 m s. They also showed that the structure of the sea breeze front is similar in 
structure to gravity currents produced in laboratory water-tanks. 

Mitsumoto et al. (1983) conducted a laboratory experiment which used a water-tank 
to simulate the sea breeze. This experiment reproduced many of the characteristics of 
the sea breeze flow, such as strong vertical motion at the leading edge of the sea breeze, 
compensating flow aloft and production of longitudinal convective vortices within the 
sea breeze flow. 

Ogawa et al. (1986) observed the gravity head structure reproduced by Mitsumoto 's 
laboratory experiment during a lake breeze. Ogawa's experiment utilized a kytoon- 
mounted ultrasonic anemometer-thermometer to measure instantaneous wind and tem- 
perature in order to observe the turbulent structure of the lake breeze. He found that 
the penetration of a lake breeze produced an upward rolling motion within the head of 
the front. He recorded upward flow velocities around 3 m s. Ogawa found that the 
temperature profile had peaks indicating the frontal passage. These peaks were due to 
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the boundary between the offshore wind aloft and the lake air below. Another feature 
that was observed was the development of the Thermal Internal Boundary Layer 
(T1BL). The T1BL develops due to the temperature discontinuity between the land and 
the water (Stunder and Sethurman, 19S5). The TIBL plays a significant role in the 
fumigation of coastal pollution that is advected into it. 

Estoque's (1962) model had a maximum penetration of sea-breeze front of 32 km 
and a minimum penetration of 18 km. The maximum penetration occurred during calm 
initial conditions and the minimum penetration occurred during offshore counter-flow. 

Pearson et al. (1983) used a model and found that with a mean flow perpendicular 
to the coast, the speed of the sea breeze front is a linear function of the onshore com- 
ponent of the mean flow. He found that the maximum vertical velocity in the front at 
any given time was not affected by the mean flow. Another characteristic that Pearson 
found was that in a negative shear case (upper flow is from land to sea), there was an 
increase in the vertical extent of the sea-breeze front. 

For a general descriptive summary of sea breeze studies, Atkinson (1981) provides 
a comprehensive review from man}' locations and experiments. 

D. LAND-SEA BRE1 ZE EXPERIMENT (LASBEX) 

A measurement program was conducted in the Monterey Bay area in September 
1987, to understand the structure of the land- sea breeze effect and generate data to 
evaluate many of the models which have been used to demonstrate the behavior of the 
air flow. This experiment was unique because it utilized a varied array of instruments 
including doppler acoustic sounder (sodar), a doppler lidar, rawinsonde launches over 
the the sea, radiosonde launches over land, conventional surface recordings of temper- 
ature, humidity and wind and corresponding observations by the NOAA 9,10 and GOES 
satellites. This combined instrument analysis allowed for comparison of the data among 
many sources. The final result will be a comprehensive view of the land-sea breeze in- 
tensity and structure. This study focuses on data from one sodar supplemented by lidar 
observations. 

The products obtained from the sodar data were the height of the marine boundary 
layer and time-height cross sections of wind direction, wind speed, vertical velocity and 
the standard deviation of the vertical velocity. From these it was possible to compute 
time-height cross sections of horizontal divergence and turbulent kinetic energy. The 
time-height cross sections will be compared with conventional data obtained from 
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radiosonde, rawinsondc and surface observations. The boundary layer height observa- 
tions will be compared with the time-height and conventional data. 

The lidar provided pictures of the radial wind field in the sea breeze front. The lidar 
observations will be compared with the sodar to find differences between the instruments 
ability to map the sea breeze front and measure the associated wind field. 

The combination of lidar and sodar has been applied before to analyzing the at- 
mosphere. However, the two instruments have not been used to study the sea breeze 
circulation. Russell et al. (1974) used the combination to observe the atmospheric 
boundary layer in an urban environment. One comparison in Russell's study was the 
ability of either instrument to detect subsidence. In the case of the lidar, detection was 
degraded because aerosols tend to be excluded from descending air. For the sodar, the 
temperature gradient associated with a large-scale subsidence inversion is easily detect- 
able. Lopez (1977) used the combination to study the properties of convective activity 
and small cumulus formation. This combination is effective since either instrument can 
detect different features to give a better overall picture of the atmosphere. 

" The objective of this thesis is to evaluate the effectiveness of the sodar in the de- 
tection and mapping of sea breeze events and to compare its measured wind field with 
that of the lidar!* A case study on 16 September 19S7 is presented. An overall look at 
typical characteristics cf the sea breeze in the Monterey Bay. Salinas Valley will be gen- 
erated for use in other studies, such as aerosol transport or the evaluation of satellite 
imagery. 
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II. DESCRIPTION OF EXPERIMENT 



A. LOCATION OF THE EXPERIMENT 

LASBEX was conducted from between 15 and 30 September 1987 on the Monterey 
Bay coast. The collection of data started at 1200 Pacific Standard Time (PST) on 15 
September and ended at 1300 PST on 30 September. All instruments were operated 
continuously except the contracted research vessel Silver Prince, which was stationed in 
Monterey Bay to monitor surface wind, temperature, humidity and pressure as well as 
launch rawinsondes. The ship collected data from 1000 PST until 1600 PST on various 
days throughout the experiment. The instruments evaluated in this paper were the 
Pennylvania State University (PSU) sodar and the NOAA Wave Propagtion Laboratory 
(WPL) lidar. The sodar was located southeast of Castroville at the mouth of the Salinas 
Valley. T he lidar was located on a flat knoll approximately 1.5 km from the coast, 
southeast of Moss Landing, it was on a line of 310° with the sodar, which approximated 
the orientation of the Salinas Valley. The separation between the sodar and the lidar 
was 5.5 km. The Silver Prince was stationed in Monterey Bay approximately 16 km from 
the lidar location. The three main instrument locations are shown in Tig. 2. 

The terrain around the sodar site was at elevations of less than 30 meters, with ele- 
vations greater than 30 meters northeast, starting 2.5 km away. The interference caused 
by the town of Castroville was considered minimal for the experiment. The terrain sur- 
rounding the lidar was generally flat with higher elevations to the north and east. To the 
east the terrain rose to levels of 100 meters and greater, up to 400 meters. To the north, 
the terrain generally had a maximum elevation of 150 meters. North-northeast of the 
lidar. Elkhorn Slough, a brackish inlet, interruped the elevated terrain. 

At the lidar, sodar and ship, twenty-second averages of vector wind velocity, air 
temperature, relative humidity and surface pressure were collected at the stations during 
the experiment. Radiosondes were launched at the sodar station at the discretion of the 
research team, coinciding with significant changes in the environment. 

Satellite AVHRR imagery during LASBEX was available from the Scripps Institu- 
tion of Oceanography. 

B. ACOUSTIC DOPPLER SODAR 

During LASBEX, the doppler sodar was used to monitor the acoustic backscattered 
energy of turbulent fluctuations in the atmosphere. The sodar is a 3-axis monostatic 
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Fig. 2. Location of Silver Prince, Lidar and Sodar: The Silver Prince was at lo« 
cation A. The Lidar was at location B. The Sodar was at location C. Ele- 



vations are in feet. Distances are in kilometers. 
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doppler system manufactured by Xontech and consists of a phased array of 25 
transducers which 'operated at a frequency of 1600 hertz. The pulse repetition interval 
is 10 seconds per beam, giving a cycle time of 30 seconds between the 3 acoustic beams. 
Two cycles were integrated per minute. A wind profile was sampled every five minutes, 
giving an effective number of 10 cycles per record. 

The sodar uses sound energy to detect the small-scale temperature fluctuations and 
discontinuities in the mean temperature profile (Crease et.al., 1977; Koracin and 
Berkowicz, 1988). The addition of doppler capabilities allows for the recording of wind 
speed and direction (Kaimal and Haugen, 1977). The backscattered energy provides in- 
formation about the height of the marine boundary layer, wind speed, wind direction, 
vertical velocity and vector winds. Neff and King (1988) demonstrated that the sodar 
can give indirect spatial and temporal information on both the location and thickness 
of turbulent la\ers. In addition, Neff was able to identify mixing processes in the 
boundary layer. Small mixing processes could be attributed to turbulence, where as 
large processes could be associated with effects such as Kevin-Helmholtz instabilities. 

C. DOPPLER LIDAR 

The lidar used during LA SB LX was the WPL 10.6 micron doppler system. The lidar 
operated under computer control and had real-time data processing and color display 
capabilites. The system could compute the zero, first and second moment of the echo 
spectrum for each range gate. For LASBEX, the lidar was operated in the range-height 
indicator (RH1) and plan postion indicator (PPI) modes. In RHI, the beam of the lidar 
sweeps in a 180° arc from foward to overhead and then behind the instrument. In PPI. 
the beam sweeps at a fixed angle of elevation from the horizon in a 360° circle around 
the lidar, the path describing an inverted cone. In both modes, the display of informa- 
tion is of the radial \elocit> relative to the lidar. Familiarity with a doppler presentation 
is required to analyze the PPI display. This analysis is similar to a doppler radar PPI 
and is described in Wood and Brown (1983). In RHI, the lidar beam was aligned parallel 
to the axis of the Salinas Valley with an orientation angle of 140°, or east-west with an 
angle of 90°, with the preferred orientation being east-west. This was because 
climatological records indicated that the sea breeze blew from due west in the Moss 
Landing area. In PPI, the elevation of the lidar beam was 1° and 5°. The 1° beam 
shows the low level structure below 200 meters, whereas the 5° beam gives structure 
above 300 meters. These values were determined triginometrically. Table 1 lists some 
of the clear-air radial velocity parameters of this instrument: 
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Table 1. W PL DOPPLER RADAR CAPABILITIES 



Expected range 


10-15 km 


Range resolution 


140-400 m 


Beam width 


< 0.1 degree 


Pointing accuracy 


0.1 degree 


Accuracy of radial velocity measurment 


> 1.0 in s 


Operation modes 


PPL VAD, RHI, range-time, raster-scan 



The lidar uses the backscattered energy of a laser from gas molecules and 
particulates to detect atmospheric features such as clouds or aerosol layers. Collis and 
Russell (1976) cited an example of how lidar was used to find the depth of the mixed 
layer or determine trie heights of imersions, which can affect the dispersion of 
pollutants. Xoonkester et al. (1974) was able to detect Kelvin-Helmholtz instabilities in 
the boundary-layer. In clear air. a lidar can detect the stratification associated with a 
deep stable boundary layer. 

D. SURFACE MEAUREMENTS 

A portable coastal meteorological station was used to take 20-second averages of 
wind direction, wind speed, air temperature, relative humidity and barometric pressure 
at the surface. This instrument used a propeller anemometer to measure wind speed. 
A self-aspirated thermometer and hygrometer measured air temperature and relative 
humidity. Wind direction was taken from a vane that gave angles with respect to true 
north. The ranges of the instrument errors were: ± 5° for wind direction: ± 0.3 m s for 
wind speed; ± 0.3° for air temperature; and ± 5% for relative humidity. 

E. SYNOPTIC SITUATION BETWEEN 16-30 SEPTEMBER 1987 

On 16 September, California was under the influence of a surface high pressure 
system off the coast and a local thermal trough in the center of the state. The winds 
were northerly to northwesterly between 5 and 10 knots. This configuration continued 
to influence the west coast from 17 to 24 September. Fig. 3shows the synoptic situation 
during this period. Around the 25th, a short wave which was detectable at the 850 mb 
level began to pass through northern California from the northwest. This wave built up 
a strong thermal gradient off the coast of California, providing an increase of cold ther- 
mal advection. This was the most significant feature associated with an actual synoptic 
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frontal passage during LASBEX. Fig. -4 shows the short wave location. Winds during 
this wave passage were mainly northerly with speeds averaging 7 m s or more. After the 
passage of the wave on the 26th, high pressure off the coast and the thermal trough in 
the center of the state were reestablished. This weather pattern persisted for the re- 
mainder of the experiment. 

The combination o ' a high pressure system over the eastern North Pacific Ocean 
and a thermal trough in central California is the typical summmer and early autumn 
synoptic pattern. This situation accounts for the common Monterey Bay pattern of low 
clouds and fog with light and variable winds during the night and morning. A tropical 
disturbance named 'Otis' formed west of the Baja Peninsula during the 23rd but did not 
significantly affect surface and lower level flow during that part of the experiment. 
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Fig. 3. Typical Surface Synoptic Situation 17-30 September 1987: Solid lines are 

sea level isobars. Values are + 1000 mb. \ 
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Fig. 4. 850mb Synoptic Situation 25 September 1987: Solid lines are isoheights. 

Values are + 1000 meters. Dashed lines are isotherms. Temperatures are 
in degrees Celsius. 
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III. DATA AND ANALYSIS 



A. SODAR THEORY 

1. Basic Sodar Theory and Turbulence in the Inertial Subrange 

This discussion of sodar theory and operation follows Underwood (19S1). The 
sodar utilizes the backscatter of acoustic energy to determine certain parameters in the 
atmosphere. The relationship between the acoustic power emitted and the backscatter- 
ing cross section is defined by the sodar equation. 



P(R) = P 0 -jf L v e~ 2 ’ R t 



\ A 0 



g( 7l), 



(3.1) 



where 



P{R) = received acoustic power (watts), 



P 0 = transmitted acoustic power (watts), 

A = cross sectional area of the receiver antenna {in 2 ), 
G = antenna gain factor. 



L v = length of the acoustic pulse in space (m). 
a = path averaged molecular attenuation coefficient (am -1 ), 
R = distance from the scatterer to the receiver (m), 



— — , - = "excess" attenuation and 

-r A t 

a{7i) = backscattering cross section (m -1 ). 

The backscattering cross section can be related to the temperature structure 
parameter O t , and the velocity structure parameter C], as shown by Little (1969). These 
structure parameters indicate the intensity of high frequency variations of temperature 
or velocity. Explicitly, the backscatter cross-section is 
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(3.2) 



o(r.) = 0.03A' ' cos*($) 



c; 

2 



2/ e 

COS- ( — 



0.13 



( sin — - ) 1 



where 

ct(tt) = backscattered acoustic energy, 
c — speed of sound, 

T = absolute air temperature, 

C] = velocity structure parameter and 
Cj = temperature structure parameter. 



Turbulent mixing in clear air causes fluctuations in potential temperature and 
water vapor which scatter and attenuate the acoustic energy in the sodar signal. The 
backscattering cross section determines how much energy is returned to the sodar an- 
tenna after attenuation by the turbulence. Turbulent eddies with a size equal to one half 
the wavelength of the emitted acoustic energy / a/2 , (Russell et al., 1974) are detectable 
with a sodar. Eddies of this size fall within the inertial subrange of turbulence and are 
assumed to be isotropic. Energy relationships in the inertial subrange follow the 
Kolmogoroff description of isotropic turbulence. The backscattered data are processed 
using East Fourier Transforms (EFT) to establish the frequency shift due to the radial 
components. The two main fields which the monostatic doppler sodar can map are the 
thermal structure and the mean velocity field along with its variances. 
a. Thermal Structure 

Molecular attenuation is described by Sutherland (1975) and Piercy et al. 
(1977) as those effects at the molecular scale which influence the propagation of energy 
through a medium such as air. The factors which cause molecular attenuation of sound 
energy are due to thermal, viscous and diffusive effects as well as rotational and 
vibrational losses on the molecular scale. The molecular attenuation increases with fre- 
quency and is the limiting factor in the maximum usable frequency in a sodar. The lower 
limit is based on antenr \ size. "Excess" attenuation is due to turbulent spreading of the 
beam pattern (Brown and Clifford, 1976). Any of the sodar signal that is not attenuated 
can be scattered back to the receiver by turbulent eddies according to (3.2). 
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